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* REMOTE  STORAGE 

Ueber  erzwungene  Wellenbewegungen  zäher, 
inkompressibler  Flüssigkeiten  in  elastischen  Röhren. 

I.  Teil. 

Historische  Einleitung. 

Die  Bestrebungen,  den  Vorgang  der  Blutzirkulation  genau 
zu  erforschen,  führten  schon  E.  Weber1)  im  Jahre  1827  auf 
den  Gedanken,  die  Pulsbewegungen  experimentell  zu  verfolgen. 
Hiebei  stellte  er  fest,  dass  diese  Pulsbewegung  im  ganzen 
Arteriensystem  nicht  momentan  erfolgt,  wie  bisher  angenommen 
wurde,  sondern  eine  gewisse,  messbare  Geschwindigkeit  besitzt. 
Er  fand,  z.  B.,  dass  die  Welle  eine  Röhre  von  35,5  mm.  Durch- 
messer, 4 mm.  Dicke  und  9620  mm.  Länge  bei  einem  Druck 
von  8 mm.  Wassersäule  in  0,964  Sek.  durchläuft,  sodass  also  die 
Geschwindigkeit  za.  10  m/Sek.  beträgt. 

W.  Weber2)  suchte  die  von  seinem  Bruder  E.  Weber 
experimentell  erhaltenen  Resultate  auf  mathematischem  Wege 
herzuleiten.  Er  bestimmte  die  Differenz  der  Flüssigkeitsvolumina, 
welche  durch  zwei  unendlich  benachbarte  Querschnitte  einer 
elastischen  Röhre,  im  Abstande  dx  voneinander,  in  dem  Zeit- 
element dt  hindurchgehen;  sie  muss  gleich  der  Vergrösserung 
des  Volumens  der  Röhre  zwischen  diesen  beiden  Querschnitten 
sein.  Ist  c die  veränderliche  Geschwindigkeit  der  Flüssigkeit, 
r der  veränderliche  Radius  der  Röhre  und  ist  c-dr<<<^r-dc  d.  h. 
die  Geschwindigkeit  sehr  klein,  so  ergibt  sich : 

de  2 dr 

dx  r dt 

Die  Elastizität  der  Röhre  wurde  dadurch  berücksichtigt,  dass  die 
Zunahme  des  Röhrenhalbmessers  der  Zunahme  des  Druckes  p 
proportional  gesetzt  wurde,  also 

dr  — adp, 


2 


£ 

wo  a = — — Vergrösserung des Röhrenhalbmessers pro lkg Druck, 

£ = „ „ „ bei  p kg.  „ . 

Diese  Annahmen  führen  auf  die  bekannte  Bewegungsgleichung 

d2  r r d2  r 

dt2  2 ae  dx2  , 

wo  q die  Dichte  der  Flüssigkeit  bedeutet. 

r 

Der  Wert  drückt  das  konstante  Verhältnis  der  beiden 

2 ao 

partiellen  Differentialquotienten  aus,  und  dieser  ist  nach  den  Ge- 
setzen der  Wellenbewegung  gleich  dem  Quadrat  der  Wellen- 
geschwindigkeit V.  Also  ist  nach  W.  Weber 


wo  M der  Elastizitätsmodul  nach  Webers  Auffassung  bedeutet. 
Derselbe  ist,  entgegen  der  gewöhnlichen  Definition,  gleich  dem- 
jenigen spec.  Druck  (Druck  dividiert  durch  die  Dichtigkeit  der 
Flüssigkeit),  welcher  nach  dem  Gesetz  der  Elastizität  einer  Ver- 
doppelung des  Röhrenhalbmessers  entspricht. 

Die  Webersche  Formel  ergibt  für  eine  Röhre  von  16,5  mm. 
Radius,  bei  einem  Druck  von  3500  mm.  Wassersäule,  eine  Ge- 
schwindigkeit von  10033  mm/Sek.  Die  unmittelbare  Messung 
ergab  eine  Geschwindigkeit  vom  11255  mm/Sek. , was  mit  der 
Rechnung  soweit  übereinstimmt,  als  bei  der  damals  erreichbaren 
Genauigkeit  erwartet  werden  kann.  Resal3)  schlägt  einen  ähn- 
lichen Weg  ein  wie  Weber  und  findet,  wenn 

v = Geschwindigkeit  der  Flüssigkeit  in  Richtung  der  Achse, 
e — Dicke  der  Röhrenwand, 

R0  = Innerer  Radius  der  Röhre, 
tc  — Dichte  der  Flüssigkeit, 
g — Beschleunigung  der  Schwerkraft, 
ds  = Abstand  zweier  unendlich  benachbarter  Querschnitte 
der  Röhre, 

E = Yongs  Elastizitätsmodul  der  Dehnung, 

d2  v Eeg  d2  v 

"dt2"  ~ 2R07t  ’ d^2"' 

Somit  folgt  für  die  Fortpflanzungsgeschwindigkeit  der  Welle 


3 


Dieser  Formel  entspricht  auch  die  Webersche,  \venn 

1 er  j*  ßg 

— — — und  — — gesetzt  wird.  Eine  genaue  Bestimmung 
q 7t  a R0 

beider  verschieden  definierten  Elastizitätsmoduln  führt  also  auf 
die  nämlichen  Resultate. 


Korteweg4)  behandelt  das  Problem  der  Berechnung  der 
Schallgeschwindigkeit  in  elastischen  Röhren.  Er  legt  die  sog. 
Scheibchenhypothese  zu  Grunde.  Dabei  wird  angenommen,  dass 
irgend  ein,  zwischen  zwei  auf  der  Röhrenachse  senkrecht  ste- 
henden  Ebenen  eingeschlossenes  Flüssigkeitsscheibchen  während 
des  Vorbeischreitens  der  Verdichtungs welle,  zwar  in  radialer 
Richtung  breiter  und  in  achsialer  Richtung  schmäler  wird,  aber 
immer  durch  Ebenen  begrenzt  bleibt,  sb  dass  von  Verbiegungen 
dieser  Ebenen  abgesehen  wird.  Ferner  wird  angenommen,  dass 
die  Wellenlänge  gross  genug  sei,  um  bei  den  in  der  Röhren- 
wand entstehenden  Spannungen  nur  auf  die  Ausdehnung  oder 
Einschnürung  des  ringförmigen  Durchschnittes,  senkrecht  zu  der 
Achse,  achten  zu  müssen,  während  die  Dehnungen  in  Richtung 
der  Achse  vernachlässigt  werden  dürfen. 

Der  von  Korteweg  eingeschlagene  Weg  ist  nun  folgender: 
Bezeichnen : 

Ri  ==  Innerer  Radius  der  Röhre, 
ai  — Dicke  der  Röhrenwand, 

Ei  = Elastizitätsmodul  der  Röhrenwand, 

E — Elastizitätsmodul  der  Flüssigkeit  (für  inkompressible 
Flüssigkeit  ist  E = ^), 

Qi  ==  Spec.  Masse  der  Röhrenwand, 

Q = „ „ » Flüssigkeit, 

x = Entfernung  eines  Flüssigkeitsteilchens  im  Gleichgewicht 
von  einer  zur  x Achse  senkrecht  stehenden  Ebene, 
w = Druck  daselbst, 
ui  = Änderung  von  x,  zur  Zeit  t, 


Pi  = 

i*i  = 

so  ergibt  sich 


w. 


Ri  „ 


n v 

Y>  n 


r, 


4 


1)  Durch  Vergleichung  des  Volumens  eines  Scheibchens  im 
Gleichgewichtszustände  und  während  der  Wellenbewegung: 


— -j- 
E 


dui 


+ 


2n 


= 0. 


dx  Ri 

2)  Durch  Berechnung  der  Beschleunigung  des  Scheibchens 
vermöge  des  Druckunterschiedes  : 

d2  m 1 dpi 

dt2  q dx 

3)  Unter  Berücksichtigung  der  Elastizität  der  Röhrenwand : 

Ei  ai  i*i 

n o P1  tV2 

d2  i*i  Ri 

dt2  ai  qi  ’ 

oder  bei  Vernachlässigung  der  geringen  Trägheit  der  Röhren- 
wand 

Ei  ai  1*1 

Pi  = 


Ri 


Unter  diesen  Voraussetzungen 
gleichung 


erhält  man  die  Bewegungs- 


d2  ui 
dx2 


+ 


2 q Ri 


d2  Ul 

dt2 


0. 


E Ei  ai 

Diese  Beziehung  gilt,  wenn  sowohl  Flüssigkeit  als  Röhrenwand 
elastisch  sind,  und  enthält  die  Fälle,  wo  entweder  nur  die  Flüssig- 
keit, oder  nur  die  Röhre  als  elastisch  zu  betrachten  ist. 

Für  inkompressible  Flüssigkeiten,  wenn  E = folgt: 
d2  ui  2 q Ri  d2  ui 
Ei  ai 


= 0 


dx2  Ei  ai  dt2 
und  hieraus  die  Fortpflanzungsgeschwindigkeit  der  Welle, 

V = 


V 


ai  Ej 


2 q Ri 

1 g 

dies  ist  wiederum  die  Resalsche  Form,  wenn  — = — gesetzt 

Q 71 

wird. 


Diese  Formel  gilt  nur  unter  den  folgenden  Bedingungen  : 

1)  Wellenlängen  gross  gegenüber  dem  Röhrendurchmesser. 

2)  Vernachlässigung  der  lebendigen  Kraft  der  Transversal- 
bewegung der  Flüssigkeit  gegen  die  der  longitudinalen. 

3)  Elastizitätsmodul  unabhängig  von  der  Grösse  der  Be- 
lastung. 
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4)  Flüssigkeit  inkompressibel,  ohne  Berücksichtigung  der 
Zähigkeit. 

5)  Vernachlässigung  der  durch  die  Biegung  der  Röhren  - 
wand  erzeugten  Längstensionen. 


Wenn  E}  den  Yong’schen  Elastizitätsmodul  der  Dehnung 
bezeichnet,  dann  ist 


Et  - Ei  1 1 — ^1  -f 


ai 


2Ri 


l -f-  f. i, 

wo  l und  fA.  die  Lame’schen  Konstanten  der  Elastizität  bedeuten. 
Setzt  man  nach  Wertheim  1 = 2 ^ so  ist 


Ei  = Ei  (l - **-). 

\ 6R,  ) 

Boussinesq5)  findet  in  seiner  Untersuchung  über  dieses 
Problem  dieselben  Resultate  wie  Korteweg. 

Eine  exakte  Behandlung  der  Gesetze  der  Schallausbreitung 
in  Flüssigkeiten,  im  Innern  elastischer,  zylindrischer  Röhren, 
ohne  Rücksicht  auf  die  Zähigkeit,  wurde  von  Lamb6)  durch- 
geführt. Die  Ableitung  der  von  ihm  gegebenen  Gleichungen  folgt 
später.  Er  findet,  dass  die  Geschwindigkeit  in  dünnen  Röhren 
zwischen 


ö 2 •, 
c0  und 


2a> 


hB 


1 4- 


2 a z 

~hB 


Co  liegt. 


Hierin  bedeuten  : 

a = Innerer  Radius  der  Röhre, 
h = Dicke  der  Röhren  wand. 

c0  — Geschwindigkeit  in  der  unbegrenzten  Flüssigkeit. 


g = — Poissonscher  Koeffizient,  wo  l und  u 

2 + |U) 

(Steifigkeit)  die  bekannten  Lame’schen  Elastizitäts- 
konstanten sind. 

E = ^ = Yong’s  Elastizitätsmodul  der  Röhre, 

X -f-  ^ 


x — Kubische  Kompressibilität  der  Flüssigkeit  = QoCo2, 
wo  = Dichte  der  Flüssigkeit. 
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Für  inkompressible  Flüssigkeiten,  d.  h.  für  sehr  grosse  z,  werden 
diese  Grenzen : 

Eh  „ Bh 

und  , 

2a^o  2a^o 

also  auch  hier,  da  sich  B nur  wenig  von  E unterscheidet,  eine 
prinzipielle  Übereinstimmung  mit  der  Resalschen  Formel. 

V.  Kries7)  gibt  ebenfalls  eine  angenäherte  math.  Theorie 
der  Schlauch  wellen.  Diese  Untersuchung  enthält  schon  die  wich- 
tigsten qualitativen  Resultate  über  den  Einfluss  der  innern  Rei- 
bung auf  die  Wellenbewegung  im  elastischen  Schlauch.  Allein 
dieser  Einfluss  ist  hier  in  ganz  allgemeiner  Weise  eingeführt, 
unter  Beibehaltung  der  Scheibchenhypothese  und  ohne  Rücksicht 
auf  die  speziellen  Eigenschaften  der  Zähigkeit  inkompressibler 
Flüssigkeit.  Soll  dieselbe  exakt  berücksichtigt  werden,  so  sind 
die  vollständigen  hydrodynamischen  Gleichungen  zu  benützen. 
Y.  Kries  geht  aus  von  der  Korteweg  entsprechenden  Formel 
(pag.  4) 

d v 1 d p 

dt  o d x’ 

ergänzt  dieselbe  durch  ein  Reibungsglied,  das  allgemein  der  Ge- 
schwindigkeit proportional  gesetzt  ist,  aber  behält  immer  noch 
die  Scheibchenhypothese  bei,  also 

d v 1 d p 

dt  G d x ^ 

Für  den  Druck  p wird  dann 

d2  p e Q d2  p o p 

dt2  o dx2  ^ dt 

Hierin  bedeuten: 

Q = Querschnitt  der  Röhre. 

x = Abstand  eines  Querschnittes  vom  beliebig  gewählten 
Anfangspunkt. 

o — Dichte  der  Flüssigkeit, 
p ===  Druck. 

q = Frequenz  in  2 u Sekunden. 

Y1 

e = Dämpfungskonstante  (für  lange  Wellen  = — also 

2 a 

sehr  klein). 
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a = Fortpflanzungsgeschwindigkeit  der  Welle. 

5-5  = «02  ist  einfach  V2  = Quadrat  der  Wellengeschwindig- 

G 

keit  nach  Resal. 

Aus  obiger  Gleichung  folgt  für  die  Geschwindigkeit  der  Welle 


1 _ J_ 

«2  ~ 2 «o2 

und  für  die  Dämpfung 


Die  Fortpflanzungsgeschwindigkeit  der  Welle  nimmt  dem- 
nach mit  zunehmender  Reibung  ab  und  ist  umso  grösser,  je 
grösser  q,  d.  h.  je  kleiner  die  Wellenlänge.  Die  Dämpfung  s 
wächst  mit  steigendem  q;  kurze  Wellen  werden  stärker  gedämpft 
als  lange. 

Für  kleine  Reibung  findet  er  ausserdem  für  die  Geschwin- 
digkeit der  Flüssigkeitsteilchen  in  Richtung  der  Röhrenachse 


Ohne  Reibung  ergab  sich 


p - A (,  x 

v ==  = e cos  q t . 


oa  oa  \ a J 

Man  sieht,  dass  (bei  sinusförmigen  Wellen)  die  Geschwin- 
digkeit der  Flüssigkeitsteilchen  mit  oder  ohne  Reibung  in 
doppelter  Weise  modifiziert  ist;  erstlich  ist  der  absolute  Betrag 
der  Geschwindigkeitsschwankungen  kleiner,  zweitens  findet  eine 
gewisse  Phasendifferenz  statt. 

Über  das  Fliessen  zäher  Flüssigkeit  in  Röhren  bestehen 
Arbeiten  von  Boussinesq8)  u.  a.  Die  nachfolgende  Untersuchung 
im  III.  und  IV.  Teil  soll  den  Einfluss  der  Zähigkeit  auf  die 
Wellenbewegung  in  elastischen  Röhren  ermitteln. 
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II.  Teil. 


Elementare  Betrachtung  der  Wellenbewegung  inkompressibler 
Flüssigkeiten  in  elastischen,  zylindrischen  Röhren 
unter  periodisch  veränderlichem  Druck,  ohne  Rücksicht 
auf  die  Zähigkeit. 


Die  Flüssigkeit  bewege  sich  im  Innern  einer  unendlich 
langen,  zylindrischen  Röhre  von  gleichförmigem  Querschnitt  und 
sei  an  der  Grenzfläche  einem  elastischen  Druck  unterworfen.  Die 
wirbelfreie  Bewegung  werde  aus  dem  Ruhezustand  durch  Kräfte 
erzeugt,  die  nur  anfänglich  wirken,  so  dass  nachher  alle  äussern 
Kräfte,  also  auch  die  Schwerkraft,  vernachlässigt  werden.  Im 
weitern  werde  angenommen,  dass  die  Bewegung  sehr  klein  sei, 
sodass  die  Quadrate  der  Geschwindigkeiten  der  einzelnen  Teilchen 
vernachlässigt  werden  können. 


Es  bezeichne : 

p = Druck  in  der  Flüssigkeit  (als  Überdruck  über  den 
im  Ruhezustand  vorhandenen  konstanten  Druck 
auf  gefasst). 

a = Elastizitätskonstante  der  Röhre  (in  gleicher  Weise 
definiert  wie  bei  Weber,  pag.  2). 

R = Innerer  Radius  der  Röhre  in  Ruhe, 
z/r  = Erweiterung  des  Röhrenhalbmessers  unter  dem 
Drucke  p zur  Zeit  t. 

q ==  Dichte  der  inkompressiblen  Flüssigkeit. 

q = Geschwindigkeit  eines  Flüssigkeitsteilchens  zur  Zeit 
t mit  den  Koordinaten  x parallel  zur  Röhrenachse 
und  dem  radialen  Abstand  r von  derselben, 
u,  w = Komponenten  dieser  Geschwindigkeit  parallel  zur 
x Richtung  bezw.  zur  radialen  Richtung. 

Wir  benützen,  da  die  Bewegung  wirbelfrei  ist,  das  Geschwin- 
digkeitspotential ; dann  ist 


d cp 

u = — — und 

B x 

wo  cp  = Geschwindigkeitspotential. 


Die  Elastizität  der  Röhre  führen  wir  ein  durch  den  ein- 
fachen Ansatz  : 


9 -r- 


(1) 


p0  -h  — J r für  r = R. 


Es  gelten  folgende  hydrodynamischen  Beziehungen  : 
1.  Die  Kontinuitätsgleichung 

(2)  n + = 


d2cp  d2(p  1 dcp 

~d72  + dr2"  ' r d7 


2.  Das  bekannte  Integral  der  Eulerschen  Bewe 
gungsgleichunge  n.9) 


(3) 


p dt  2 


wo  F(t)  eine  willkürliche  Funktion  der  Zeit  darstellt,  die  in 


dcp 

dt 


eingeschlossen  werden  darf. 

Für  p setzen  wir  den  Wert  aus  Gl.  (1)  in  Gl.  (3)  ein,  die 

konstante  Grösse  — wird  ebenfalls  in  ^ eingeschlossen,  dann 

q dt 

folgt,  da  q2  verschwindend  klein,  und  wenn  — =x  gesetzt  wird: 

a q 

(4)  Jr  = — — für  r = R. 

x dt 

Bei  den  verschwindend  kleinen  radialen  Amplituden  der 
Bewegung,  wird  die  Normalkomponente  der  Geschwindigkeit  der 
Flüssigkeitsteilchen  der  Normalkomponente  der  Grenzfläche  in 
erster  Annäherung  gleich  sein,  somit : 

dqp  6(J  r) 

dr  dt 

Differentieren  wir  Gl.  (4)  nach  t und  addieren  Gl.  (5),  so  folgt 

(6) 


(5) 


w = 


für  r = R. 


+ * £2  = 0 für  r = R. 
dt2  dr 


Für  stehende  Wellenbewegungen  mit  der  Frequenz  - ° = n und 

2 TC 

2 TC 

der  Wellenlänge  l = setzen  wir 

(7)  <p  = P cos kxe1  ((Tt  + f) 

wo  P nur  eine  Funktion  von  r darstellt. 

Diesen  Wert  für  cp  führen  wir  in  Gl.  (6)  ein  und  bekommen 
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(8) 


_ 8cP 


(J-1  (p  =z 


dr 


für  r = R. 


Aus  der  Kontinuitätsgleichung  (2)  erhalten  wir  durch  Einsetzen 
von  (7) 

d2  p i q p 

-) P k2,  oder  weil  — k2  — (ik)2 


(9) 


dr2 

d2P 


r dr 
ikr 


dP 


+ P = 0. 


d (ikr)2  ikr  d (ikr) 

Dies  ist  die  bekannte  Form  der  Besselschen  Differen- 
tialgleichung.10) Hieraus  ergibt  sich  als  allgemeines  Integral 

P = A J (ikr)  -f  B Y(ikr), 

o « 

wo  J (ikr)  = Besselsche  Funktion  nullter  Ordnung  1.  Art  und 

Y (ikr)  = » » » » 2.  » 

In  der  Röhre  darf  P natürlich  nie  unendlich  gross  werden. 
0 ° 

co,  so  muss  das  Glied  mit  Y (ikr)  ver- 


Da  aber  lim.  Y (ikr) 
r — 0 

schwinden,  d.  h.  es  muss  die  Konstante  B 
deshalb  als  Integral  der  Gl.  (9) 


o sein.  Es  bleibt 


(10)  P = A J (ikr). 

Diesen  Wert  in  Gl.  (7)  eingesetzt  ergibt: 

(11)  <P  — A J(ikr)  cos  kxe1  (<rt  + f)  . 

Unter  Berücksichtigung  dieser  Gleichung  folgt  nach  Gl.  (8) 

(12)  -2_  ..  kiJ(ikR) 


o = — X 


J (ik  R) 


Dieser  Ausdruck  ist  reell,  da  i J (ik  R)  reell  ist. 

Nach  Gl.  (4)  ergibt  sich  für  die  Erweiterung  des  Röhren- 
halbmessers 


(13) 


J r = i 


. AJ(ikR)-u  _ n ^ ((jt  + 6) 


cos  (kx)  • eM 


oder,  weil  für  uns  nur  der  reelle  Teil  in  Betracht  kommt, 
(13a)  z/  r = 


J (ik  R)  cos  (kx)  sin  (ut  -f-  e). 


x 


Diese  Gleichung  stellt  ein  System  stehender,  radialer  Wellen 
dar.  Um  ein  System  fortschreitender  Wellen  zu  erhalten,  su- 
perponieren  wir,  analog  dem  Vorgang  in  Lamb9)  für  Oberflächen- 
wellen, 2 Systeme  stehender  Wellen  von  der  gleichen  Wellen- 
länge. Dabei  müssen  sie  einen  Phasenunterschied  von  7*  Wellen- 
länge aufweisen,  damit  die  Berge  und  Täler  des  einen  Systems 
mit  den  Knoten  des  andern  zusammenfallen.  Dann  erhalten  wir 
auf  analoge  Weise  wie  oben 

(14)  <p=  A J(ikr)ei(ört-kx) 

und 

(15)  Jr=  i — j(ikR)eil,Tt-k’t) 

X 

oder  als  reeller  Teil 

(15a)  Jv  = — — J(ikR)sin(ot-kI). 

% 

2 7t  2 7t 

Es  ist  die  Wellenlänge,  T — — die  Schwingungsdauer, 

k g 

somit  folgt  für  die  Fortpflanzungsgeschwindigkeit  der  Welle  aus 

l — c T,  c = — , also  nach  Gl.  (12) 
k 


V 


— i J (ik  R) 


\ 


i J(ikR) 


k J(ikR)  \ a^k  J(ikR) 

Dieser  Ausdruck  gilt  für  beliebige  Wellenlängen.  Wird  aber  X 
im  Verhältnis  zu  R sehr  gross,  dann  wird 

— i J (ikR)  = Und  j (ik  R)  — 1}  also 


(17) 


/ z R / R 

Vk^l  .'“V  * 2 


in  Ueberein- 


R . /JL 

l v " 2 V 2a^’ 

Stimmung  mit  Weber  (pag.  2),  also  auch  mit  Resal,  Korteweg,  etc. 
Unsere  Berechnung  zeigt  somit,  wie  weit  dieser  Ausdruck  bei 
kleiner  werdenden  Wellenlängen  zu  modifizieren  ist. 
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III.  Teil. 

Wellenbewegung  zäher  inkompressibler  Flüssigkeiten 
in  dünnen,  zylindrischen,  elastischen  Röhren. 

§ 1.  Die  hydrodynamischen  Gleichungen. 

Die  allgemeinen  Gleichungen  für  zähe  Flüssigkeit,  ohne 
Einwirkung  äusserer  Kräfte  lauten  nach  Lamb9): 

'du  , du  , du  du 

■ -f-  w — 
dz 

/dv  dv  dv  dv' 

\ dt  dx  dy  dz 

dw  . dw\ 


(18) 


e - 


U\  A 5P 

-NM»-/ 

7 / d X 


= ^Av-|£ 

dp 
dz 


/dw  . dw  , - ~ ..  , 

<Mt:+u7-  + vi_  + wt  =-“Aw  - 

\dt  dx  dy  dz/ 

Hierin  bezeichnen  u,  v,  w die  Geschwindigkeitskomponenten 
parallel  zu  den  Koordinaten  Achsen  in  einem  Punkte  (x  y z)  zur 
Zeit  t.  Für  irgend  einen  bestimmten  Zeitmoment  t geben  sie 
die  Bewegung  in  irgend  einem  Punkte  des  von  Flüssigkeit  er- 
füllten Raumes,  wo  p — Druck,  q = Dichte  der  Flüssigkeit, 

d2  d2  d2 

.f.i  = Reibungskoeffizient  ist  und  Adas  Symbol  für 1 

dx2  dy2  dz2 

bedeutet. 

Da  die  Geschwindigkeit  sehr  klein  vorausgesetzt  ist,  können 


~ . du  du 

die  Grossen  u — , v — , 
dx  d y 


du  . ...  du 

w — etc.  gegenüber  — , 
dz  dt 


d v 

— etc.  ver- 
dt 


nachlässigt  werden,  und  es  bleibt 

du  A dp 
Q — = p Au—  — 


dt 


dx 


(19) 


dv  dp 

dt  dy 

dw  A dp 

e — = mAw—  — 

dt  dz 


Da  die  Bewegungserscheinungen  symmetrisch  zur  Röhren- 
achse sind,  so  führen  wir  die  Zylinderkoordinaten  x und  r ein,  wo 
die  x-Achse  mit  der  Röhrenachse  zusammenfällt.  Bezeichnet 
jetzt  u die  longitudinale  Geschwindigkeitskomponente  längs  der 
x-Richtung  und  w die  radiale  Geschwindigkeitskomponente,  so 
wird  nach  Stokes  J1) 


Wir  führen  die  Stokes’ sehe  Strom funktion  ip  ein, 
indem  wir  setzen 


1 3 ih  1 d ib 

— und  w = 4 

r d r r d x 


(21)  u 
Damit  wird  die  Kontinuitätsgleichung 

(22)  — -f-  — w -|“  — - — 0 erfüllt. 

3 x r dr 

Setzen  wir  in  Gl.  (20)  — = v (wo  v der  sog.  kinematische 

Reibungskoeffizient  ist),  eliminieren  den  Druck  p,  indem  wir  die 
obere  Gl.  nach  dr  und  die  untere  nach  ox  ableiten  und  die 
obere  von  der  untern  subtrahieren,  so  ergibt  sich  in  symbolischer 
Schreibweise 
_d_ 

dt 


(23) 

Nach  (2 


dx2  Br2 


wird 

d w du 
<3x  dr 

also  wird  Gl.  (28) 

d / d2 


1 /d2ip  d2ip  1 B ip\ 
r l dx2  dr2  r dr  ' 


dt 


t 77,  + 7;  + i-- 


oder 

(24) 


[t 


dx2  dr2  r dr 


B2ip  d2\p  1 B ip 


ld 

v 

[_dt 


dx2  dr2 
d2  d2 


r dr 

1 d 


H 


dx2 


d r2  r d r 


' d2 
dx2 


d2 


dr2  r dr  r2 


^)=0. 


Zur  Abkürzung  sei  symbolisch  — (-  — — — — = D 

oxj  d r2  r 3 r rj 

gesetzt,  dann  lautet  Gl.  (24) 
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(25) 


V I) 

et 


KH 


Setzt  man  nach  Stokes  xp  = t/q  -f-  t/q,  wobei  D = 0, 

ist 

(£-’d)(dtH^dI(£-'d)t  -a 


SO 


Die  Gleichung  (25)  wird  also  auch  erfüllt,  wenn 


a)  D 


xpi_ 

r 


= 0 und  b)  I 


\et  J r 


0. 


Für  ipi  und  rp2  bestehen  somit  die  folgenden  Differential- 
gleichungen : 

^20)  a‘2^i  i &V*  1 


(27) 


dx 2 

e2xp2 


Sr2 

S2lp2 


r <9  r 
1 


0 

l e i{j2 


e X" 


5 r 


r e 


et 


0. 


Wir  suchen  ein  partikuläres  Integral  dieser  Gleichungen, 
das  die  Erscheinungen  einer  einfachen  Sinuswellenbewegung  mit 

vorgeschriebener  Frequenz  n =-^—  darstellt,  setzen  also  sowohl 

2 7t 

ip1  wie  auch  t/q  proportional  e1  + kx)  . Es  sei  t/q  = P . el  (crfc  + kx)  , 
wo  P nur  eine  Funktion  von  r,  dann  wird  Gl.  (26) 


(28) 


— — k3P  = 0 oder,  weil  — k2  = (ik)2 , 

d r2  r d v 

52P  1 aP  (ik)2  P = 0. 


dr2'  r <9r 
Hieraus  folgt 12) 

P = r [Ai  J(ikr)  -f-  EYx(ikr)]  also 

weil  auch  hier  lim.  Y1  (ikr)  = oo,  muss  die  Konstante  E,  aus  den- 
r ==  o 

selben  Gründen  wie  früher,  Null  sein.  Es  wird  daher 
t/q  — r A:  J1  (ikr)  ei((7t  + kx) 

Aber  da  i J1  (ikr)  reell  ist,  setzen  wir  Ai  = iA,  also 
(29)  t/q  = rAiJ1(ikr)ei(ört  + kx) 

Da  nach  obigem  auch  t/q  = P'el((rt  + kx),  wird  Gl.  (27), 

— k2  = (ik)2  gesetzt  wird, 


wenn 
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= o 


d2P' 


l dp' 


r dr 


(ik)2P' 


1 G 


P' 


oder,  wenn 

(30) 

(31) 


— + k2  = /?s 

v 


gesetzt  wird, 


d2P'  1 ap' 


+ (i^p'. 


d r2  r d r 

Diese  Form  entspricht  genau  derjenigen  der  Gl.  (28),  deshalb 
folgt 

P'  = r B J1  (ißr) 

also,  wenn  wir  statt  B,  analog  wie  bei  (29),  i C setzen, 

(32)  rCi  J1(i/Jr)ei<fft+k*> 

Da  ip  = i pl  ip2,  so  wird 

(33)  ip  — v [Ai  Jx(ikr)  -f-  CiJ^ißr)]  eUo-t  + kx) 
wo  ß durch  Gl.  (30)  gegeben  ist. 

Durch  diesen  Ausdruck  für  ip  wird  die  Differentialgleichung 
(25)  erfüllt.  Gl.  (33)  ist  also  ein  Integral  derselben. 

Um  eine  Beziehung  zwischen  der  Stromfunktion  ip  und  dem 
Druck  p zu  finden,  gehen  wir  aus  von  Gl.  (20),  in  welcher  wir 

— — v und  für  u und  w die  Werte  aus  (21)  einsetzen;  dann  er- 
gibt sich: 

(34) 


i e2ip 

V 

-( 

' d2  ip 

d2ip 

, 1 d^ 

1 1 £lf. 

r dr  dt 

r 

dr\ 

K d x2 

8 r3  1 

r 8 r ) 

q d x 

i e2ip 

V 

A_  / 

/ o2  . 

, 8 *P 

1 dtp'' 

i _I£p 

r dxdt 

r 

öx' 

v 1 dx2 

dr2 

r 8 r , 

' q d r 

Durch  Addition  der  Gl.  (26)  und  (27)  erhält  man 

d2ip 


(35). 


d ip  .1  d p 1 d ip2 

' r ör  v dt 


d x2  d r2 

Diesen  Wert  setzen  wir  in  Gl.  (31)  ein  und  bekommen 


(36) 


1 d2ifj 

v d / 

1 8rp2\ 

l _I£P 

r drdt 

1 

> >■  8 1 j 

' Q d X 

1 d2  i p 

_v__8_( 

, l5^ 

v dt  1 

l_  lap. 

r dxdt 

r 8 x l 

' e dr 

’i  + so  folgt 
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(37) 


^p=e  oV  und  Sv  = _s 

dx  T di  dt  6 r r d x dt 


Diese  Beziehungen  finden  sich  bereits  in  der  erwähnten 
Arbeit  von  Stokes11). 

Es  sollen  p und  \p  proportional  el(ört  + kx)  sein.  Dadurch 
beschränken  wir  unsere  Untersuchung  auf  rein  stationäre,  fort- 
schreitende Wellen  mit  bestimmt  vorgeschriebenem  Schwingungs- 
zustand für  p. 

Multipliziert  man  beide  Gl.  (37)  mit  r,  übt  auf  die  Erstere 

c)  d 

- - — und  auf  die  Zweite  aus,  addiert  beide,  so  folgt: 

0 x or 


(38) 


o p 


o2 

o p 


1tSi=0 


(39) 


+ p = O. 


d x2  d r2  1 r d r 
Für  p = P"(r)ei((7t  + kx)  gesetzt,  folgt  aus  (38)  wenn  — k: 
= (ik)2  wieder  die  bekannte  BesseFsche  Differentialgleichung 
d2 P"  1 d P" 
d(ikr)2  ikr  d (ikr) 

Somit 

P"  = CJ(ikr)  + DY(ikr). 

Auch  hier  muss  D — - 0 sein,  damit  p nicht  für  r = o un- 
endlich wird;  es  ist  daher: 

(40)  p = CJ(ikr)ei(at  + kx) 

Unter  Berücksichtigung  von  Gl.  (37)  u.  (29)  wird 

(41)  p = — A ga  J (ikr)  ei(at  + kx)  . 

Führt  man  in  den  Gl.  (21)  für  ip  den  Wert  aus  Gl.  (33) 
ein,  so  erhält  man  für  die  Geschwindigkeitskomponenten  u und 
w die  folgenden  Ausdrücke : 

A kJ  (ikr)  -f-  Cß J (i/?r) j e1 

- Ak  J1  (ikr)  - Ck  JHi/*r)j  e‘  (nt 

p = — A q o J (ikr)  e1  (ff  1 + kx) 

= r^AiJ(ikr)  + CiJ(i/?rj  | ei(<It  + “ 


u = 


w 


wo 


(a  t + kx) 


+ kx) 
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ß=±  \/k2+v- 

Da  aber  k2  für  lange  Wellen  (wie  sie  für  unsere  Unter- 
suchung allein  in  Betracht  fallen)  sehr  klein  und  v klein,  also 

— gross  gegenüber  k2  , so  kann  man  setzen 
v 

p = ± V1  7 = ± Vv ^ qVf 

o i i k r 

Ferner  kann  für  k2  <I<C  1,  J(ikr)  = 1 und  J(ikr)  = gesetzt 

2 

werden,  dann  folgt : 


II. 


(ort  -f  kx) 


WO 


u = j^Ak  -)-Cq\/iJ  (rq  \J — i)j  e1 
w = -^Ak2r+  C kJ1  (rq  \/—~  ij|  e1 


p = - A Qae' (<rt  + 

Tp 


= r j — + C i J1  (rq  \J—  i)  j e1  <<r  ‘ 


(a  t + kx) 


+ kx) 


§ 2.  Die  Gleichungen  der  Elastizität. 


Wir  greifen  ein  Volumelement  der  Röhren  wand  heraus,  be- 
zeichnen es  mit  J V und  setzen  voraus,  dass  es  so  dünn  sei, 

2 
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dass  die  Variation  der  Wanddicke  vernachlässigt  werden  könne. 


Bezeichnet : 

R = Innerer  Radius  der  Röhre  im  Ruhezustand. 

J cp  = Winkel,  unter  dem  das  Element  von  der  Röhrenachse 
aus  erscheint. 

h = unveränderliche  Dicke  der  Röhrenwand. 

J x — Länge  des  herausgegriffenen  Elementes. 

u'  = Verschiebung  in  der  Längsrichtung  bei  Einwirkung 
eines  radial  gerichteten  Druckes  p. 

w'  = Verschiebung  in  radialer  Richtung  unter  demselben 
Druck  p. 

•p  = Spannung  in  Richtung  der  x-Achse. 

£ — Spannung  in  Richtung  des  Umfanges. 

E = Yongs  Elastizitätsmodul  der  Röhrensubstanz. 

O — Poissonscher  Koeffizient  ==  Koeffizient  der  Querkon- 
traktion zur  Längendilatation.  (Diese  von  Lamb  ab- 
weichende Bezeichnung  benützen  wir,  um  die  Ver- 
wechslung mit  der  Frequenz  o zu  vermeiden.) 

Q0  ==  Dichte  der  Röhrensubstanz. 

Es  ist  z/V  ===  r - Jtp  hJx,  somit  die  Masse  dieses  Elementes 
z/M  = q0  • r • z/ cp  • hz/x. 


1 


— und  — Q sind  die  oben  genannten  Werte  für  die  Dicke  1, 
h h 

so  dass  sich  ergibt : 


Dehnung  in  der  Längsrichtung  : 


Dehnung  in  tangentialer  Richtung : 
Hieraus  folgt : 


(42) 


und 
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(43) 


j^e/v  au: 

1 — <92\R  d x 


z/  x ) h • R ^99  — • h • R • Jcp  = J x • R • z/  cp. 


Wirkt  aber  in  der  Längsrichtung  in  irgend  einem  Quer- 
schnitt die  Spannung  ^3,  so  erreicht  diese  Spannung  im  Abstande 

d vn 

z/x  von  diesem  Querschnitt  den  Wert  $ -] J x.  Die  Kraft- 

d x 

Wirkung  auf  unser  Volumelement  in  der  Längsrichtung  beträgt 
demnach 

. asp/h 

h dx 

Somit,  da  die  Beschleunigung  in  longitudinaler  Richtung 
d2u'  . „ , , , d2u 

ist,  Q0  • R • J cp  • h • J x = — d x • R J cp  also 

dt2  dt2  dx 

(44)  eo  h — — = — . 

dt“  d x 

Da  der  innere  Überdruck,  welcher  in  Richtung  wachsender  r 
wirkt,  p ist,  so  wirkt  auf  das  Röhren element  der  Druck  p • R • 
J(p-Jx;  in  entgegengesetzter  Richtung  wirkt  infolge  der  Ela- 
stizität die  Kraft—  Jx-h  mit  ihrer  Komponente  — Öz/xh-  sin  — ■ 
h h 2 

— Q-Jx-hJp.  Somit,  da  die  Beschleunigung  in  tangentialer 
Richtung  ^ W 


dt2 

Q0  • h • R • J cp  J x 
(45) 


d2w' 


= p-R-z/</)-Jx — £}-Jx-h*z/^,  also 


Po 


d2w' 

"dt2" 


— p D 

r 


Substituiert  man  Gl.  (42)  in  (44)  und  (43)  in  (45),  so  erhält 
man  die  Lambschen  Gleichungen.*) 


(46) 


wo  B = 


d2 


u 


B /d2  u'  , Öd  w'  \ 


dt2 

d2w' 


Qo 

P 


dx2 


E 


dt2  h^0 


Qo 


R d x 
G du' 


und 


w 


R dx  R2 


1 — G2 


*)  Genau  genommen  müssen  noch  die  durch  die  Zähigkeit  der  Flüssig- 
keit erzeugten  Druckkomponenten  mitberücksichtigt  werden.  Statt  p wäre 
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(47) 


Da  u',  w'  und  p proportional  e1((7t  + kx),  so  ergibt  sich  : 

2 B \ , i 6 B , 

2 u -j — w = 0 

q0  k R q0 

1 B\  # p 


und 


O i B 
— • — u'  + 
kR  <>0 


w 


k2  R2  Q0  J k2hQo 

In  den  hydrodynamischen  Gleichungen  treten  die  Geschwin- 
digkeitskomponenten auf  und  nicht  die  Verschiebungen  wie  in 


Gl.  (47);  wir  setzen  also 


du ' 


u und 


w 


==  w,  folglich  sind 


6 1 5 t 

diese  Gleichungen  noch  nach  der  Zeit  t zu  differentieren.  Sie 
lauten  dann  : 

6 » B 

R 


(48) 


W = 0 


und 


Qo 


0i 

kR 


B_ 

^0 


U + C2  + 


1 


k2R2 


w = 


iap 
k *hQo 


§ 3.  Bed ingungsgleichungen  für  die  Geschwindigkeit 

langer  Wellen. 


Die  Gleichungen  II  charakterisieren  die  Bewegung  der 
Flüssigkeit  für  lange  Wellen.  Die  Gleichungen  (48)  dagegen 
diejenige  der  elastischen  Röhre.  An  der  Röhrenwand,  für  r = R, 
müssen  die  Geschwindigkeiten  u und  w der  Gl.  II  denjenigen 
der  Gl.  (48)  entsprechen.  Wir  führen  die  Werte  aus  II  in  (48) 
ein,  sondern  die  Konstanten  A und  C ab,  dann  folgt  aus  (48) 

A • a -j-  C ß — o 
A y Cd  = o, 


(49) 


p — 2 /li zu  setzen,  ebenso  ist  in  der  tangentialen  Richtung  ein  Zusatz- 

5 r 

glied  zu  D mit  dem  Faktor  ^ beizufügen.  Diese  Glieder  können  aber  in 

g 

den  Gl.  (46)  im  Verhältnis  zu  p und  zum  grossen  Faktor  — vernachläs- 

Qo 

sigt  werden,  um  so  mehr  als  eine  ev.  Reibung  in  der  elastischen  Wand 
sowieso  nicht  berücksichtigt  wird. 
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wo 


(50) 


« = Rk(D  -j-  — 


ß = J(Rq\/—  i)  | R q V/i  D ^ )L _Jl_  j Qh  ] 

l .T(RnV/_il  J 


J (Rq\/— i) 


y = i 0 b i 


k2  R2  cs 


. b . . ru2R  , . . 

l hi (wo  ö=  k • c) 

2 ‘ k2  V 


a = J(Rq\/—  i)|  q \/i  — ©b  y 
l k j 


J (Rq\/ — i)  * 


f Rk  (V 


J(Rq  \!  i) 

b 

k2  R2 


B 


E 


wenn 


= b,  c2  — b = D und 


= % gesetzt 


e e(i—  ©3)  ' e»h 

wird.  Da  die  Determinante  der  Gl.  (49)  null  sein  muss,  so  er- 
halten wir  als  Bedingungsgleichung  für  k 
(51)  a • ö = ß • y. 

Ferner  ergibt  sich  aus  Gl.  (49)  der  Wert  für  die  Kon- 
stante C 

A • a 


(52) 


C = — 


ß 


Führt  man  in  Gl.  (51)  die  Werte  (50)  ein,  berücksichtigt, 

dass  q\/i  =4~q\/ — h und  setzt  R 9 \! — i ==  Z,  so  ergibt  sich 
i 

o 

J(Z) 


- DZ  j 

1 J(Z) 


b / k2  R2 

i0b^i<9b-i^-  + ic2('r  + 


(58) 

= |^.Z^+Rk(b-^- 

kR  j(Zb  V VW 


Rk(c2-  b + |b 


Aber  r = — ist  bei  nicht  zu  dicken  Röhrenwänden  relativ 
Qo  h 

k2  R2 


gross,  also  kann  für  lange  Wellen 
lässigt  werden. 


2 


gegenüber  t vernach- 
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J (Zl 

Weil  y- — - eine  komplexe  Grösse  ist,  setzen  wir 


J(Z) 


(54) 


(55) 


Rq  v/—  j(RW^i)  = z m = P + iQ. 


J (Rq  V—  i) 

Setzen  wir  ferner  zur  Abkürzung 

1 


J(Z) 


0 


S = b - 

v2 

so  ergibt  Gl.  (53) 

(56)  [D(P  + iQ)  + 0bj  [k2S 


L = 2 D + 0b, 


Rl7(72] 

+ k2[(P  + iQ)0b 


b] 


L 


0. 


Dies  ist  die  Bedingungsgleichung  zur  Berechnung  von  k. 
Da  k eine  komplexe  Grösse  ist,  setzen  wir  k = m -f-  in 

wo  m — / = c T — Wellenlänge,  T = Schwingungs- 

X (7 

dauer,  n = Dämpfungskonstante  in  Richtung  der  Fortpflanzungs- 
geschwindigkeit. 

Den  Wert  für  k2  — (m2  — n2)  -f-  i 2 mn  und  o — m • c in 
Gl.  (56)  eingeführt,  gibt : 

— rm2  c2[R D (P  + iQ)  -j-  R 0b]  + (m2  - n2  + i2mn) 
(57)  .[D(P  + iQ)  + 0b]  S 

-(-  (m2  — n2  -f-  i 2 mn)  [(P  i Q)  0 b — • b]  — ==  0. 


m 


Diese  Gleichung  gestattet  zunächst,  das  Verhältnis  — zube- 

n 

rechnen.  Zerlegen  wir  sie  in  ihre  reelle  und  imaginäre  Kompo- 
ponente,  ordnen  nach  den  Variabein  n und  m,  so  folgt: 
m2  a'  — n2  ß'  — 2 mn/'  = o 
m2  d'  — n2  / -f-  2 mn/?'  = o, 

wenn  zur  Abkürzung  folgende  Setzungen  gemacht  werden  : 


D 


Pb  \ P 

— - 7C2RP-  b)  + 02b2  I 1 

2 / V 2 


z c2  R 0 b 


ß'  =-(p 

2 


2)  (D  + 02  b) 


/ = Q 
d’  = Q 


(D  + & 


2 b)  | 


D / - — T c2  R 
2 


©2b2 


2 , 


(59) 
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Dividiert  man  die  Gl.  durch  mn,  addiert  beide  und  ordnet, 
so  folgt : 

m ß'2  _L_  a/2 

— = 2- — — : die  obigen  Werte  eingesetzt,  ergibt: 

n a'r'  — ß'd' 

rn  = 2 b (c2  — b 02  b)  (Q2  + (P  - 2)2 ) *) 
n 2Q[2(c2— b)  + 0b]  rc2R 

Wir  berechnen  ferner  aus  Gl.  (56)  die  Fortpflanzungs- 
geschwindigkeit c.  Zu  dem  Zweck  schreiben  wir  sie  wie  folgt : 


(61) 


|D(P  + iQ)  + öbJ ^2S-2rR^ 
= — j^(P  -)-  i Q)  0 b — b j L. 


Wir  setzen  links  den  Ausdruck 

(62)  2 r R — = X + iY. 

k2 

Die  reellen  und  imaginären  Teile  der  Gl.  (61)  sind  dann  : 

1.  Reeller  Teil : 

(63)  — X [D  P + 0b]  + Y D Q = 

= — 2S[DP  + 0b]  - (P  0b  — b)L 

2.  Imaginärer  Teil : 

(64)  XDQ  + Y[D-P+  0 b]  = Q0bL  -f  2 S • D • Q. 

Aus  diesen  beiden  Gleichungen  folgt  durch  Elimination  von  Y 

X = |"(DP  + @b)  { — 2 S(D  • P + 0b)  — (P  0b  — b)  • L } 


(65) 


(66) 


— DQ2{©bL  + 2S-D 


(DP+  0b)2  + D2  Q2 
Löst  man  die  Klammern  auf  und  ordnet,  so  wird 

x _ b (D  + 02  b)  [(P  - 2)  (D  • P + 0b)  + Q2  D] 
(D  • P + 0b)2  + D2  Q2 


ist 


*)  Anmerkung:  In  den  für  uns  in  Betracht  kommenden  Fällen 

m o 

— sehr  gross.  Denn  für  <?  = = 1,  b eo  1 0 , 9 — 0,5,  Rr  = 10, 

c2  106  (s.  pag.  35 — 36)  und  den  Werten  für  Rq  = 10,  P cv»  Q cs;  10  (s. 

pag.  34)  ergibt  sich  — -1?  • ^ rvj  102  . 

n 10  . 108 . 106 . 10 

Es  darf  also  jedenfalls  n2  gegen  m2  , in  erster  Annäherung  aber  auch 
n gegen  m vernachlässigt  werden. 
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X ist  der  reelle  Teil  von  2rR  — . Für  k = m-(-in  wird,  wenn 

k2 

n2  <C<<  m2  (’s.  Anmerkg.  pag.  21),  X = 2 t R c2. 

Dies  in  Gl.  (66)  eingesetzt,  ergibt  in  erster  Annähe- 
rung, wenn  c2  zunächst  gegen  b vernachlässigt  wird,  also 


B 


Qo 


und  D + S2  b = — b (1 


e2) = — — 


h-E  (P  — 2) (P  - ©)  + Q: 

L 


c2K. 


D = — b 
gesetzt  wird  : 

(67) 

2R  Q (P  — Of  + Q2  _ 

Ist  die  Reibung  v = o,  also  q = 4L  = oo,  so  wird  auch 

P = Q — oo,  also  K = 1 und  es  bleibt 
2 h • E 


(68) 


Co 


2 R q 


, die  von  Resal  gefundene  Form  (s.  p.  3). 


Der  Ausdruck  K 


(P  — 2)  (P  - 0)  + Q' 


b e r ü cksich- 


(P  - 0)2  + Q2 

tigt  die  Zähigkeit  der  Flüssigkeit,  er  ist  also  der 
Resalschen  F ormel  als  Korrektionsfaktor  beizufügen. 

Für  eine  genauere  Berechnung  kann  der  so  korrigierte 
Wert  in  Gl.  (66)  eingeführt  werden.  Es  wird  dann  in  2.  An- 
näherung : 


D — c — b 
H = 1 


2 TT 

c0  K 


b = — b H 


wo 


4k-1 


h (i —& 
2 R 


K dann  wird 


(69) 


C = C0 


2H-©2  (P  — 2)  (HP  — ©)  + H Q2  2 


Co  Ko, 


1 — 03  (HP  — 0)2  + H2Q2 
wo  K0  ein  verbesserter  Korrektionsfaktor  ist.  Wenn  nötig,  kann 
durch  fortgesetzte  Substitution  in  Gl.  (66)  der  Wert  von  c2 
mit  beliebiger  Genauigkeit  ermittelt  werden. 


§ 4.  Die  Longitudinale  und  die  radiale 
Geschwindigkeit  der  F 1 ü s s igk  e i t s t e il  c h e n 
bei  einfach  harmonischer  Schwingung  von  p. 

In  Gl.  (52)  ist  die  Konstante  C durch  die  Konstante  A aus- 
gedrückt. Da  C komplex  ist,  so  schreiben  wir  zweckmässig,  da- 
mit die  Ausdrücke  für  die  Geschwindigkeiten  u und  w möglichst 
einfach  ausfallen, 
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(70) 


C = 


A k (<£  + i !F)  Ak  (<2>  4-  i V)  (—  \/ - i) 


q\/i  q 

wo  0 und  lF  zu  bestimmende  Funktionen  von  R sind.  Aus  Gl. 
(52)  folgt,  wenn  für  a und  ß die  Werte  (50),  sowie  die  Substitu- 
tionen L und  Z Gl.  (54)  und  (55)  eingeführt  werden. 

(71)  & + = LRc^A 


2J1(Z)[DRq\/ 


- J(Z) 
JHZ) 


i Ob] 


Nach  unserer  Setzung  (54)  wird  im  Nenner 

o o 


Rq  \! 


r-J(Z) 


J(Z) 


R 

i — =■=  — q V --  l 

JHZ)  i 4 v JHZ) 


(P  i Q)* 


Führen  wir  nun  die  Thomsonsche  12)  Bezeichnung  ein 
(72)  — \/ — i JL(Z)  — ber'  Z i bei'  Z,  so  wird 

# -f-  i = — Rq\/r  -pL = 


V'-1  2[ber'Z  + i 
iRq 


ibei'Z]  | - 


Hr+iQ)-i©b 


oder 

(73)  0 + i = 

2 [ber'  Z + i bei'  Z]  D(P  + iQ)  + Ob 

Setzen  wir  noch  zur  Abkürzung 

M DP  -|-  Ob  und  D Q = N, 

so  erhält  man  aus  (73)  durch  Trennung  des  reellen  und  imagi- 
nären Teiles 

Rq  L (M  bei'  Z -f-  N ber'  Z) 


(74) 


0 = 


ip  = _Rq 


[ber'  Z2  bei'  Z2  ] (M3  + N2) 
L (M  ber'  Z — N bei'  Z) 


2 [ber'  Z2  + bei'  Z2  ] (M3  + N2 ) 

Wird  in  erster  Annäherung  c2  gegenüber  b vernachlässigt, 
also,  weil  D = c2  — b,  L = b (0  — 2),  M = b (O  — P)  und 
N = — b Q gesetzt,  so  folgt 

Rq (2  - 0)[(P  — O)  • bei'Z  + Q ber'Z] 


(74a) 


0=  — 


— 


2 [ber'  Z2  + bei'  Z2].[(P  — O)2  + Q2) 
Rq  (2  — 0)[(P  — 0) ber'Z  — Qbei'Z 


2 [ber'Z2  + bei'  Z2]-[(P  — O)2  + Q2] 
Mit  der  Substitution  (70)  ergeben  die  Gl.  II : 
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Erstens:  für  die  Geschwindigkeit  u in  Richtung  der 

Röhrenachse,  wenn  wir  rq\/ — i — z setzen  und  nach  Thomson 

o , _ 

(75)  J (r  q y — i)  — her  z -|-  i bei  z einführen, 

(76)  u = Ak[l  — (0  -f  i W)  { berz  + i bei  z }]ei(<rt  + kx) 

Da  k = m -(-  i n,  wo  aber  n nach  Gl.  (60)  gegenüber  m 
vernachlässigt  werden  darf  (s.  Anmerkung  auf  pag.  23)  folgt, 

— nx 

(77)  u = Ae  m [1  — (<Z>  + i ¥) (berz  + ibeiz)] e1  ((T  1 + mx) 

Da  wir  im  Folgenden  nur  die  Schwingungen  der  Flüssig- 
keitsteilchen betrachten,  die  in  der  Gleichgewichtslage  in  einem 
bestimmten  Querschnitt  (z.  B.  x = o)  liegen,  so  braucht  auf  diese 
Dämpfung  nicht  mehr  weiter  Rücksicht  genommen  zu  werden, 

— nx 

so  dass  wir  den  Faktor  e weglassen.  Er  ergibt  sich  als  reeller 
Teil  aus  Gl.  (77) : 

(78)  u = A m [(1  -[-  <£> berz — lFbeiz)cos/  — (*fberz-| -$beiz)sin/], 
wo  / — g t -f-  m x 

Bezeichnen  wir  den  Koeffizient  von  sin/  mit  E und  den- 
jenigen von  cos  yv  mit  F,  so  ist 

u = — Esin/  -f-  F cos/  — Gcos(/  -j-  <5),  also  G = -f-\/E2  -f-  F2 
somit : 


(79) 


Am 


V(1 


<2>berz  — iFbeiz)2-|-(^berz+  <£beiz)2-cos(/-|-d), 


wo 


E 


tg  d = — -- 

F 1 


*Fberz  4-  beiz 


-f-  & ber  z — lP  bei  z 
Zweitens:  für  die  Geschwindigkeit  w in  radialer  Rich- 
tung, ebenfalls  unter  Vernachlässigung  der  Dämpfung  n, 


(80)  w 


-H«  + i «0(— V— Ü JHrq \/— i)]  e,(ot 


+ mx) 


und,  wenn 

(81)  ( — \/ — i)  JL(rq  \/ — i)  = ber'z  + i bei'z, 
ergibt  der  reelle  Teil  von  (80) : 

Am2  T 

(82)  wT  = ((Dber'z — W bei'z)  cos/ 

q I 


-f W ber'z  -f-  & bei'  z)  sin  / 
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oder,  analog  behandelt  wie  u : 
Am2 


(83) 


w 


wo 


^ -|-  ^ber'z  -f~  ® bei'zj  -|- 


+ j (D  ber'z  — ^Fbei'zJ  cos  (x  -f-  ül), 


tgd' 


qr 

~2 


-f-  llJ  ber'z  -|-  ( D bei'  z 


0)  ber'z  — <Fbei'  z 


§ 5.  Untersuchung  der  Bewegung  bei  Annahme 
rein  radialer  Schwingung  der  Röhren  wand. 

Die  bisherigen  Ausdrücke  zeigen,  dass  stets  u,  die  Geschwin- 
digkeit in  Richtung  der  Röhrenachse,  ausserordentlich  viel  grösser 
ist  als  die  Geschwindigkeit  w senkrecht  zur  Röhrenachse.  Die 
Röhrenwand  macht  also  hauptsächlich  longitudinale  Schwin- 
gungen. Betrachten  wir  deshalb  den  speziellen  Fall,  wo  sie  nur 
radiale  Schwingungen  vollzieht. 

Es  sei  also  in  der  Röhre  für  r = R die  Geschwindigkeit 
der  Flüssigkeit  an  der  Röhren  wand  in  Richtung  der  Röhren- 

dw 

achse  u = o;  ebenso  setzen  wir  — . da  dieser  Betrag  sehr  klein 

dx 

ist,  gleich  null. 

Aus  Gl.  (46)  folgt  unter  dieser  Annahme  und  unter  Be- 
rücksichtigung der  Bemerkung,  die  zu  der  GL  (48)  führte,  sowie 

B n 

= b, 


mit  der  Substitution 
d2  w 


Qo 


1 dp 


w (für  die  Röhren  wand). 


(84) 

dt2  h^o  dt  R2 
p und  w setzen  wir  auch  hier  proportional  el((7t  + kx),  dann  ist 

(85)  — (T2  w 
Es  ist  <j2  gegenüber  sehr  klein,  wir  lassen  es  also  weg  und 

Rj2 

bekommen  als  Bedingungsgleichung  für  k : 


1 . b 

i ap w. 

ho0  F R2 
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(86) 


b . o *) 

— w = i p. 

R2  h^0  F 


Da  wir  an  der  Röhrenwand  r = R,  die  Geschwindigkeit  u = o 
annehmen,  folgt  ans  II 

,_o 

A k -J-  C q \i  J (Z)  = o,  hieraus 

Ak 


(87) 


C = 


q\/i  J (Z) 

Diesen  Wert  für  die  Konstante  C führen  wir  im  Ausdruck  für 
w in  II  ein  und  bekommen  für  r = R unter  Berücksichtigung 
von  (54) 

(88)  /lR  lR 


w = — A k2  I — — 


2 P + iQ 

Diesen  Wert  in  (86)  eingesetzt,  gibt 


ai  (ff  t + kx) 


(89) 


k2 


1 


g2  R 


\ 2 P -j-  i Q, 

Aber  k2  = (m2  — n2)  -J-  i 2 m n,  o2  = m2  c2  und 

1 P . Q 


P + iQ  P2  + Q2  P2  + Q2 

Also  wird  Gl.  (89): 

(90)  (m2 — n2  -)-  i 2 mn)  ^ ^ i : Q 


+ i 


2 P2  + Q2  P2  + Q 
n R 


hieraus  ergeben  sich  für  die  reelle  und  imaginäre  Komponente 
folgende  Gleichungen : 


(91) 


(m2  — n2)  ( — — ^ 


2 P2  + 


— 2 mn 

2 2 ^ 

— nr  + t — = o 


Q 


P2  + Q: 


(m2  — n2)  — — 1-  2 mn  ^ ^ 


P2  + Qs 


2 P2  + Qs 


= o. 


*)  Bemerkung:  Nach  Gl.  II  wird  für  v=o,  wenn  keine  Reibung 

vorhanden  wäre,  w — — — A k2  R e1  ((y  t + kx)  und  p = — Aq  oel  + kx) 

Da  a — k . c,  folgt  aus  Gl.  (86) 

c2  = — --0 ; dies  ist  wieder  die  Resalsche  Form,  wenn 
2Ro 

B E 

b = — = — gesetzt  wird,  was  in  erster  Annäherung  statthaft  ist. 

Qo  (?• 
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Aus  diesen  beiden  Gleichungen  berechnet  sich  das  Verhältnis 

1 


m 

— zu: 
n 


(92) 


m 


= — 2 


(P2  -f-  Q2  — 2 P)2  + Q< 
(P2  + Q2)Qc2rIT~ 


Eine  Überschlagsrechnung  ergibt  der  Grössenordnung  nach  den 
nämlichen  Wert  wie  Gl.  (60),  so  dass  wir  auch  hier  zur  Berech- 
nung von  c,  n gegenüber  m vernachlässigen  können.  Unter  dieser 
Voraussetzung  ergibt  die  1.  Gl.  (91) 


2 1 (P2  + Q- 

nr  — 


2P) 


|2  r»2 


R 


nrc^r  — ==  o 


und 


P2  + Q2 

hieraus,  wenn  durch  m2  dividiert  wird,  und  die  Werte  von  b und 
x (pag.  20)  eingesetzt  werden, 


(93)  c2=  — •-(P2+Q2~ 
cR  2 (P2  + Q- 


2P)  Eh  (P2  + Q2  — 2P) 


2R^  (1 

Eh 


G2)(P2  + Q2) 


2Ro 


K', 


wo  K'  das  Korrektionsglied  angibt,  das  der  Resalschen  Formel 
beigefügt  wird. 

Zur  Berechnung  der  Geschwindigkeit  u bei  radialen 
Schwingungen  der  Röhrenwand  setzen  wir  den  Wert  der  Kon- 
stanten C aus  (87)  in  II  ein  und  bekommen : 


(94) 


u = Ak(l  — — — ) e‘  ‘ + kx>. 
V J(Z) 


Für  k = m -|-  in  und  mit  der  Thomsonschen  Bezeichnung  nach 
(75)  und  unter  Vernachlässigung  der  Dämpfung  in  der  x-Rich- 
tung,  wird  der  reelle  Teil  in  (94),  wenn  ot  -|-  mx  = %, 

berZ-berz-f-  beiZ  beiz' 


(95) 


u = Am 


berZ2  + beiZ2 
berZ-beiz  — beiZ-berz 


cos  / -f- 


berZ2  -f-  beiZ2 

oder  analog  behandelt  wie  Gl.  (79) 


sm  % 


]• 
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(06) 


Am-cosfe  -f-  d) 
berZ2  -f-  beiZ2 


^ [ber  Z- beiz — bei  Z • berz]2  -\- 


wo 

tg<3  = 


-f-[berZ2+  beiZ2  — {berZ-berz  -f-  bei Z- beiz}]2, 
ber  Z • beiz  — bei  Z • ber  z 


berZ2  — beiZ2  — [berZ  • berz  -]-  beiZ  • beiz] 

Auf  gleiche  Weise  ergeben  sich  aus  II,  nach  Einführung  der 
Gl.  (87),  für  die  Geschwindigkeit  in  radialer  Richtung 
die  folgendest  Gleichungen: 

r J*(z)  1 


(97) 


w — — A k2 


Bi 


e*  (O  t + kx) 


2 q\/iJ(Z)J 

Unter  den  nämlichen  Bedingungen  wie  früher  wird  die  reelle 
Komponente 


(98)  w = Am2R 


oder 


T ber 

L (t 


Zber'z  4-  beiZ  bei'z 


(berZ2  -f-  beiZ2)Rq' 
berZbei'z  — beiZber'z 


cos  x 


— | 

(berZ2+ beiZa)Rq  2RJ 


sin  x 


(99) 


Am-  • cos(x  -f-  ö ) / 

w = i / 

q(berZ2  -f*  beiZ2)  \/ 

ber  Z bei'z  — bei  Z • ber'z  — 

— ^berZ2  -f-  beiZ2) 

H 

j^berZ ber'z  -J-  beiZ  bei'z  j 

WO 


berZbei'z  — beiZber'z  — (berZ2  beiZ2) 


berZber'z  -f-  beiZbei'z 

§ 6.  Untersuchung  der  Bahn  eines  einzelnen 
Flüssigkeitsteilchens. 

Es  seien  (x,  r)  die  Gleichgewichtskoordinaten  eines  Flüssig- 
keitselementes in  einer  Meridianebene,  (xr)  die  Koordinaten  der 
Verschiebung  aus  dieser  Mittellage.  Nun  ist: 

^chp 

r d r 


(100) 


dx d (x  -f-  cc) _ 

dt  dt 


d r d (r  -f  r)  . 1 S xp 

— - = — — w = d — -. 

dt  dt  r dx 
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Für  u und  w setzen  wir  die  gefundenen  Werte  aus  den  Glei- 
chungen (79),  (83),  (96)  und  (99)  ein  und  bezeichnen  zur  Ab- 
kürzung die  Amplituden  von  u und  w mit  f bezw.  F,  dann  ist 


dx  = f cos  (ut  -f-  m (x  -j-  x)  -f-  ö)  dt 
dr  = F cos  (<rt  -f~  m (x  -f-  x)  -f-  d')  dt. 


(101)  j 


Bei  der  Integration  nach  t können  im  Argument  die  mit 


dem  kleinen  Faktor  m multiziplierten  Werte  x als  konstant 
angesehen  werden,  also 


x — f — sin  (et  m x -f  dj 
r = F — sin  (<rt  -f-  m x -f-  d'). 


(102) 


Aus  diesen  Gleichungen  ergibt  sich  nach  Elimination  der 


Zeit  t: 

(103)  x2  F2  -f-  r2  f — 2 xr  ■ Ff  cos  (d'  — d)  — ~ f2  F2  sin2  (d' — ö). 


Dies  ist  die  Mittelpunktsgleichung  einer  Ellipse  mit 
variabler  Neigung  der  Achsen. 

Zur  Bestimmung  dieser  Neigung  gegen  die  Röhrenachse 
transformieren  wir  die  Gl.  (103)  auf  ein  neues  Koordinatensystem, 
indem  wir  setzen: 


x = x'  cos  y — r'  sin  y 
r =?  x'  sin  y -j-  r'  cos  y,  woraus  folgt 


xf2  [F2  cos2  y -{-  f2  sin2  y — 2 F f cos  (d' — d)  cos  y sin  y] 


-f-  r'2  [F2  sin2  y -{-  f2  cos2  y 2 F f cos  (d' — d)  cos  y • sin  y] 


- f2  cos  y sin  y -j-  F f cos  (d' — d) 
F2  sin2  (d' — d)  = 0. 


Um  die  Normalform  herzustellen,  muss  die  letzte  Klammer 
Null  sein,  also: 


(106) 
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Da  die  Amplituden  der  Geschwindigkeit  w senkrecht  zur 
Röhrenachse  stets  sehr  viel  kleiner  sind  als  diejenigen  von  u in 
Richtung  der  Röhrenachse,  also  F « f,  so  ist  aus  GL  (106) 
ohne  weiteres  ersichtlich,  dass  die  grossen  Achsen  der  Ellipsen 
der  einzelnen  Flüssigkeitsteilchen  praktisch  parallel  zur  Röhren- 
achse liegen. 


IV.  Teil. 

Diskussion  der  erzwungenen  Wellenbewegung  einer  zähen 
inkompressiblen  Flüssigkeit. 

§ 1.  Zusammenfassung. 

In  unseren  Untersuchungen  des  III.  Teiles  stellten  wir  die 
hydrodynamischen  Gleichungen  für  die  Bewegung  einer  zähen, 
inkompressiblen  Flüssigkeit  auf  und  suchten  unter  Benützung 
der  Stokes’schen  Stromfunktion  ein  Integral  für  Wellenbewe- 
gungen im  Innern  einer  unendlich  langen,  zylindrischen  Röhre 
mit  kreisförmigem  Querschnitt,  ohne  Rücksicht  auf  die  Schwer- 
kraft. Hiebei  ergab  sich  Gl.  (33). 

Die  Grenzbedingungen  der  Flüssigkeitsbewegung  folgten 
aus  der  Elastizitätstheorie  mit  der  Annahme  einer  elastischen 
Röhrensubstanz  und  unter  der  Voraussetzung,  dass  die  Röhre 
im  Vergleich  zum  Durchmesser  nur  eine  geringe  Wandstärke 
besitze.  Dabei  wurde  jedoch  die  Reibung  im  Innern  der  Röhren- 
wand unberücksichtigt  gelassen,  obwohl  sie  möglicherweise  einen 
merklichen  Einfluss  auf  die  Fortpflanzung  der  Wellenbewegung 
der  Flüssigkeit  ausüben  könnte. 

Die  ganze  Betrachtung  beschränkte  sich  auf  einen  statio- 
nären Zustand;  auf  die  Entstehung  der  Wellen  wurde  nicht 
Rücksicht  genommen.  Die  erzwungenen,  fortschreitenden  Wellen 
haben  eine  bestimmte  gegebene  Frequenz  o (in  2 7z  Sek.)  und 

haben  eine  sehr  grosse  Wellenlänge  X — so  dass  also 

k2  « 1 ist.  Für  den  periodisch  oszillierenden  Druck  wurde  eine 
einfache  cos.  Schwingung  angenommen,  p = — A q g cos  (j  t k x), 
wo  q die  Dichte  der  Flüssigkeit  bezeichnet. 
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Die  Dämpfung  in  Richtung  der  Röhren achse  ergab  sich 
aus  der  komplexen  Grösse  k = m -f-  i wo  n die  Dämpfungs- 
konstante bedeutet.  Wie  sich  aus  der  Ueberschlagsrechnung 

Seite  21  zeigte,  ist  das  Verhältnis  — von  einer  zu  vernach- 
° rn 

lässigenden  Grössenordnung. 

Unter  diesen  Bedingungen  gelten  unsere  Gleichungen,  die 
hier  zusammengestellt  sind. 


I.  Fortpflanzungsgeschwindigkeit  der  Wellenbewegung. 

a)  Für  allgemeine  Wellenbewegung. 

In  erster  Annäherung: 


= f)  + Q2  = c0  VK, 

V 2 R e v (P  — &f  + Q2  0 v 


Gl.  (68) 


w 


o cQ  = ^ die  bekannte  Formel  von  Resal  ist,  oder 


JE 
R Q 

in  zweiter  Annäherung: 


/h  • E /H—  G2 

C “ V 2 R P V 1 — e2 


C * D 

wo  H = 1 K und  b = — = 


— 02  (P  — 2)  (IIP  — ©)+HQ2 


(HP  — 0)3  + H3Q2 


= C„V/K 


05 


E 


Gl.  (69) 


b q q (1— O2) 

b)  Für  rein  radiale  Schwingung  der  Röhrenwand. 


/ hE  /P2  + Q‘‘ 

V 2 R Q V (1  — ©2] 


— 2 P 


J)  (P'  + Qf) 


= c0  ^K' 


Gl.  (93) 


II.  Geschwindigkeiten  der  Flüssigkeitsteilchen. 

a)  Für  allgemeine  Wellenbewegung. 

1)  In  Richtung  der  Röhrenachse: 

u = Am  ^1  (1-J-^berz — ¥ bei  z)2-f-  (lFberz  -f-  0 beiz)2  cos  (%-f  d), 

Gl.  (79) 


ber  z -(-  0 bei  z 
1 -f-  ber  z — lF bei  z* 


3 
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2)  In  radialer  Richtung  senkrecht  zur  Röhrenachser 
Am2 


~_\/  S£-f-  *Fber'  z -f-  $bei'  zj  ber'z — *Fbei' zj 

• cos  {x  + <J)r 


wo 


tg  d'  = 


^ ber'  z -j-  0 bei'  z 

0 ber'  z — bei'  z 


wo 


Q = — 


¥=  — 


Rq  (2  — 6>)[(P  — (9)  -bei'Z  -f  Qber'Z] 
2 * [ber'Z2  + bei'Z2]  [(P  — Sf  -f  Q2]' 

Rq  (2  — Q)  [(P  — 9)  ber'Z  — Q bei'Z] 

2 [ber'Z2  + bei'Z2]  [(P  — 0)2  + Q2]‘ 


Gl.  (83> 


Gl.  (74a). 


u 


b)  Für  rein  radiale  Schwingung  der  Röhrenwand. 

1)  In  Richtung  der  Röhrenachse: 

_ Am  cos  (x  -f-  <5) 


berZ2  + beiZ2 


- <5)  / 

~ w [ber Z- beiz  — bei Z- ber z]"  -f- 


-f-  [berZ2  beiZ2  — { berZ  • berz  -f-  bei  Z • beiz  }]2 
wo  Gl.  (96); 

^ ber  Z- beiz — beiZ -berz 

g ~ ~ berZ2  — beiZ2  — [berZ  • berz  -f  beiZ  • beiz]' 


2)  In  radialer  Richtung,  senkrecht  zur  Röhrenachse 
Am2  cos  (%  <?') 


w 


wo 


q(berZ2  + beiZ2) 


v/ 


berZ  • bei'z  — beiZ  • ber'z  — 


( berZ2  -f  beiZa)]2+  jberZ  • ber'z  + beiZ  • bei 


2 V 


‘”1 


Gl.  (99). 


berZ  • bei'z  — beiZ  • ber'z  — (berZ2  -\-  beiZ2) 
ber  Z • ber'z  -}-  bei  Z • bei'z 


t»i'  = 
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III.  Die  Elongationen  x,  r der  Flüssigkeitsteilchen. 

1 


x = f — sin  (ö- 1 — j—  mx  -f  d) 
o 

r = F — sin  (a  t + mx  4 d') 
0 


tg2  y = 


2Ff  -cos(d'  — d) 


Gl.  (102) 


Gl.  (106) 


F2  — f2 

wo  f = Amplitude  bezw.  Maximalwert  von  u 
F ==  Amplitude  bezw.  Maximalwert  von  w 
y =.  Neigung  der  Achse  der  Bahnellipse  zur  Röhrenachse. 

Ferner  bedeuten  in  den  obigen  Ausdrücken 
A = Eine  Konstante,  die  sich  aus  p — A q g cos  [a  1 4"  kx)  be- 
stimmen lässt. 

E = Yong's  Elastizitätsmodul  der  Dehnung. 

0 = Poissonscher  Koeffizient  (Verhältnis  der  Querkontraktion 
zur  Längsdilatation). 

Q0  = Dichte  der  Röhrensubstanz, 
h ==  Dicke  der  Röhrenwand. 

R = Radius  der  Röhre. 
q — Dichte  der  Flüssigkeit. 

q = — wo  v — kinematischer  Reibungskoeffizient  (Visko- 

sität) der  Flüssigkeit,  o — Frequenz  in  2 jt  Sekunden. 
Z = Rq\/ — i und  z = rq\ — i. 

° 

J (z)  = ber  z 4"  i hei  z = Besselsche  Funktion  nullter  Ordnung 

I.  Art. 

— \j — i J1  (z)  = ber'z  4~  ibei'z  ==  Besselsche  Funktion  I.  Ord- 
nung I.  Art  nach  Thomsons  Darstellungen. 

J(Z)  V.  !) 


P4-iQ  = Z 


JHZ) 


Rq\/- 


J'(Rq\/-i) 


=4,  (Gl.  54.) 


§ 2.  Berechnung  der  Funktionen  P und  Q, 
die  durch  Besselsche  Funktionen  dargestellt 
sind,  und  Zusammenstellung  der  Konstanten. 
Wir  benützen  nach  obigem  die  von  Thomson  13)  eingeführte 
Bezeichnung  ber,  für  den  reellen  und  b e i für  den  imaginären 
Teil  der  Besselschen  Funktion  : 
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J (rq  \/ — i)  = ber  (rq)  — (—  i bei  (rq) 

— \J — i J1  (rq  \/ — i)  = ber'(rq)  i bei  (rq), 


wo 


ber(rq)  = 1 — 


(rq) 

2242  1 22426282 


+ 


(rqf 


(rq) 


12 


22426282102122 


+ 


bei  (rq)  == 


(rq)6  , (rq); 

224262  ^ 22426382102 


v10 


Es  sei 


Rq  Rq 


2^ — K und  RqV-i  = 2 x \J2  \/ — i = x ; 

ferner  ist  nach  Gray  and  Mathews14) 


xj(x)_  x («/a  f m 

2 J4(x)  2 12 

11  (x/2) 


CM 

48 


(x/2)8  13(x/2)1 


180 


8640 


12 


-(‘+4+-) 


26880 

Deshalb  folgt  näherungsweise  für  z 1 

Rq  \/ — i J (R q \J—  i) 

2 Jx(Rq  -1) 

+ i*2  (l 

Für  % }>  1 gibt  Sommerfeld15)  näherungsweise  den  Wert 

Rq  \J — i J (Rq  \J — i)  / , _1  , 3 

2 _ JHRqy/^  ~ \ I_h64z 

»_ 

V 64  z 128  z2 

Für  10  > z > 1,5  gilt  sehr  genau  nach  Zenneck16) 

RqV~1  J°(Rq^~*l  = (0,997  x + 0,277)  + i (1,007  x — 0,040). 
2 J (Rq  y— i) 

Multiplizieren  wir  diese  Gl.  mit  2,  so  folgt  für  unsere  Gl.  (54) 
Rq  ]/■ 


i •yRqL=^  = P + iQ=2  (0,997  x + 0,277) 
J (Rq  V— i) 


+ i 2 (1,007  x 
Also  wird: 


0,040)  für  1,5  < * (=^§  ) < 10- 
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P = 2fl  + Ö+-) 


P — 2 
P = 2 


3-64 
q R 1 

W + " 

0,997 


4^64qR/ 


qR +0,277 


und 


Q_  qaRa/  q‘R4 


fqR 

Q = 2L7t 


3y/8 
64  qR 


3 • 8 


128  q2R2 


für  z «<  1 
für  z 1 
für  10  > z > 1,5 

für  z < 1 
für  z 1 


•- 


Q 


1,007 

2|yrqR 


0,04oj  ft 


für  10  +>  z +>  1,5 


In  der  folgenden,  Jahnke  und  Emde  12)  entnommenen  Tabelle, 
sind  für  verschiedene  Argumente  Rq  die  numerischen  Werte  von 
P und  Q,*)  sowie  von  her,  bei,  ber'  und  bei'  angegeben. 


Rq 

p 

Q 

ber  (Rq) 

bei  (Rq) 

ber'  (Rq) 

bei'  (Rq) 

0 

2,0000 

0 

1,000 

0 

0 

0 

0,5 

2,0000 

0,0626 

0,999 

0,063 

- 0,008 

0,249 

1,0 

2,0002 

0,2494 

0,984 

0,249 

— 0,062 

0,499 

1,5 

2,0516 

0,5554 

0,921 

0,558 

- 0,210 

0,730 

2,0 

2,1610 

0,9612 

0,752 

0,972 

— 0,493 

0,917 

2,5 

2,3494 

1,4272 

0,399 

1,457 

— 0,944 

0,998 

3,0 

2,6360 

1,9016 

- 0,221 

1,938 

- 1,569 

0,881 

3,5 

2,9840 

2,344 

— 1,194 

2,283 

— 2,336 

0,435 

4,0 

3,3556 

2,746 

— 2,563 

2,293 

— 3,135 

— 0,491 

4,5 

3,7256 

3,116 

— 4,299 

1,686 

— 3,754 

— 2,053 

5,0 

4,0860 

3,474 

— 6,230 

0,116 

- 3,844 

- 4,354 

5,5 

4,4380 

3,832 

— 7,974 

— 2,790 

— 2,907 

- 7,373 

6,0 

4,7874 

4,186 

— 8,858 

— 7,335 

— 0,293 

— 10,846 

8,0 

6,1912 

5,628 

+ 20,974 

- 35,017 

+ 38,294 

- 7,662 

10,0 

7,5880 

6,874 

+ 138,841 

+ 56,370 

+ 51,373 

+135,230 

*)  Der  Hälfte  dieser  Werte  P und  Q entsprechen  in  der  Tabelle  von 

w'  Li'  w 

Jahnke  und  Emde  Seite  147,  die  Ausdrücke  — bezw. 

wo  Wo 
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Ueber  die  Grösse  der  in  Flüssigkeiten  und  Röhrensubstanzen 
auftretenden  physikalischen  Konstanten  geben  die  folgenden 
Tabellen  Aufschluss. 


1.  Viskosität  und  Dichte  einiger  Flüssigkeiten.17) 


Flüssigkeit 

Viskosität 

= 

beider  Temp. 
t° 

Dichte 

Q 

Wasser 

0,0178  -r-  0,0028 

0 — 100 

1 

Terpentinöl 

0,0146  -P  0,007 

20  — 80 

0,87 

Olivenöl 

0,808  — 0,115 

20  — 80 

0,91 

Glyzerin 

42  — 7,8 

2,8  — 20,9 

1,24 

Blut*) 

0,07  — 0,108 

37° 

1,053-1,066 

2.  Elastizitäts koeffizienten.17) 


Substanz 

Elastizitätsmodul 
in  Dyn/cm2 

E 

Poissonscher 

Koeffizient 

0 

Dichte 

Qo 

Eisen 

21000  • 108 

0,243-0,310 

7,8 

Kupfer 

10500  • 108 

0,348 

8,9 

Blei 

1700  • 108 

0,375 

11,4 

Glas 

4700—7900  • 108 

0,22-0,31 

2, 4-3, 8 

Holz 

500-1000  • 108 

— 

0,5 

Kautschuk 

0,02—0,8  • 108 

0,5 

0,92—0,99 

Für  die  numerischen  Beispiele  wählen  wir  die  Konstanten 
so,  wie  sie  etwa  einer  Röhre  aus  Kautschuk  entsprechen,  in 
welcher  eine  relativ  zähe  Flüssigkeit,  wie  z.  B.  Blut,  mit  einer 
Frequenz,  die  ungefähr  der  Pulsfrequenz  entspricht,  Wellenbe- 
wegungen ausführt. 


Es  sei  q — Q o 


= i)b=A 

Qo 

0,38  • 108. 


E 


Q (izr^)=0’5,1°8DWcm2 
(7  möge  zwischen  1 und  10, 


0 — 0,5  also  E 
v zwischen  0,01  und  0,1,  R zwischen  1 und  10  cm  variieren. 


Da  diese  3 Grössen  immer  in  der  Form  Z = Rq  = R 
zusammen  auf  treten,  so  können  natürlich  für  bestimmte  Werte 

*)  Nach  E.  Münzer  und  F.  Bloch18)  ist  die  Viskosität  des  normalen 
Blutes  ca.  4 — 6 mal  grösser  als  für  Wasser. 
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von  Z noch  viel  mannigfachere  Kombinationen  von  R,  er,  v vor- 
genommen werden. 

Das  Verhältnis  der  Dicke  der  Röhrenwand  zum  innern 
Radius,  h/R,  sei  im  allgemeinen  0,1  also  Rr  = 10. 

§ 3.  Fortpflanzungsgeschwindigkeit  der  Wellen  und 
Berechnung  der  Konstanten  und  lF. 

Unter  den  angenommenen  Verhältnissen  ergibt  sich,  ohne 
Rücksicht  auf  die  Zähigkeit,  nach  Resal: 

c» = V tu!  = 1369  cm/sek- 

Der  Korrektionsfaktor  erster  Annäherung,  der  den  Einfluss  der 

Zähigkeit  ergibt,  wird  für  Z = R = 10,  also  z.  B.  für 

eine  Röhre  von  1 cm  Radius,  bei  einer  Frequenz  a — 10  (in 
2 7t  Sek.)  und  einer  Zähigkeit  v — 0,1,  K = 0.891,  somit 
c = c0  \j  K = 1290  cm/sek.  Es  ist  ersichtlich,  dass  die  Zähig- 
keit einen  nicht  zu  vernachlässigenden  Einfluss  auf  die 
Fortpflanzungsgeschwindigkeit  der  Welle  ausübt.  Die  Formel 
für  die  erste  Annäherung  ist  in  unserem  Fall  genügend  genau, 
obwohl  hier  c2  (oo  1,7  • 106)  gegenüber  b (oo  0,5  • 10b)  ver- 
nachlässigt wurde;  denn  der  Korrektionsfaktor  der  2.  Annäherung 
wird  K0  — 0,885,  woraus  sich  c zu  1285  cm/sek  ergibt,  d.  h. 
ein  Wert,  der  von  dem  Vorigen  nicht  zu  unterscheiden  ist.  Die 
Art  und  Weise,  wie  c mit  der  Zähigkeit  variiert,  zeigt  die  nach- 
folgende Tabelle  und  Figur  (2)  für  den  Korrektionsfaktor  K, 
wobei  auch  der  Korrektionsfaktor  K'  für  rein  radiale  Schwingung 
der  Röhrenwand  angegeben  ist.  Beide  Kurven  nehmen  zu,  wenn 
Z wächst,  wenigstens  innerhalb  der  hier  angegebenen  Werte. 
Die  Fortpflanzungsgeschwindigkeit  der  Welle  nimmt  also  mit 
abnehmendem  Röhrendurchmesser,  abnehmender  Fre- 
quenz und  mit  zunehmender  Viskosität  ab.  Der  Ein- 
fluss des  Erstem  macht  sich  aber  viel  stärker  geltend,  da  er 

gegenüber  im  Quadrat  auftritt;  der  Einfluss  der  Frequenz  und 

der  Zähigkeit  dagegen  macht  sich  in  gleich  starker  Weise  gerade 
entgegengesetzt  geltend.  Ferner  folgt,  dass  in  erster  Annäherung 
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die  Elastizität  der  Rölirenwand  keinen  wesent- 
lichen Einfluss  auf  die  Fortpflanzungsgeschwin- 
digkeit ausübt;  denn  in  den  beiden  Korrektionsfaktoren 
K u.  K'  tritt  nur  der  Poissonsche  Koeffizient  auf,  welcher  sich 
für  die  verschiedensten  Substanzen  nur  in  verhältnismässig  ge- 
ringen Grenzen  bewegt,  während  b nicht  auftritt. 


R<y=  01  23456789  10 


z= 

Rq 

K 

K' 

0 

0 

0 

0,5 

0,0018 

0,0013 

1,0 

0,027 

0,020 

1,5 

0,143 

0,140 

2,0 

0,324 

0,305 

2,5 

0,493 

0,503 

3,0 

0,608 

0,675 

3,5 

0,680 

0,782 

4,0 

0,728 

0,859 

4,5 

0,760 

0,915 

5,0 

0,786 

0,958 

5,5 

0,805 

0,990 

6,0 

0,821 

1,018 

8,0 

0,865 

M oo 

10,0 

0,891 

1,141 

Vergleichen  wir  noch  im  folgenden  die  Formel  von  v.  Kries 
(pag.  7)  mit  unsern  Resultaten. 

Die  Formel  von  v.  Kries  lautet  in  unserer  Bezeichnungsweise 
du  1 d p 


d t 


Q 


— 7]  U. 


Wir  gingen  aus  von  den  Eulerschen  Gleichungen 

du  ldp  /d2u 


d2  u 


22 

Ö ü+  

B y2  + S z2 


etc.  etc. 
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Nehmen  wir  an,  dass  u nicht  von  r abhängig  sei,  legen 
also  die  Scheibchenhypothese  zu  Grunde,  so  ist 


d_u 
d t 


1 d p 
Q d x 


r>2 

O U 


+ v . 2> 


dau~e1[5,  + kI\-:  = -kJu,  folgt 


’ o 2 

dx 
du. 
d t 


1 d p 2 

— — V k u. 

Q O X 


Diese  Gleichung  ist  mit  derjenigen  von  v.  Kries  identisch,  wenn 
ri  —k2v;  es  ist  also  rt  selber  eine  Funktion  der  Frequenz. 

O m 2 V O 2 

Setze  k = m = ~,  dann  ist  y\  = v m = — hierin  kann 

c 

für  c die  Resalsche  Wellengeschwindigkeit  c0  gesetzt  werden,  so  dass 
v o 2 

7]  — — 2 (wo  Co2  = «o2  von  v.  Kries) 
co 

Dies  in  die  v.  Kries’sche  Formel  für  die  Fortpflanzungsgeschwin- 
digkeit eingesetzt,  ergibt: 


also 


c 


(für  Scheibchenhypothese) 


Der  v.  Kries’sche  Korrektionsfaktor,  der  sich  schreiben  lässt 

1 2 2 

1 £ £ ° , zeigt  auch  eine  Abhängigkeit  von  v,  von  g und 

4 co2 


v/ 


von  R,  abei*  doch  in  wesentlich  anderer  Weise,  als  in  unseren 
Formeln.  Der  Einfluss  der  Zähigkeit  der  Flüssigkeit  ist  also 
durch  den  v.  Kries’schen  Ausdruck  nicht  richtig  wiedergegeben. 
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Die  Konstanten  0 und  lF  sind 
gegeben  durch  die  Gl.  (74a)  wenn  c2  gegen 
b zu  vernachlässigen  ist.  Die  Werte  von 
0 u.  lF  finden  sich  in  nebenstehender  Ta- 
belle ; der  graphische  Verlauf  für  ver- 
schiedene Rq,  also  für  verschiedene  Röh- 
rendurchmesser, verschiedene  Frequenzen 
und  verschiedene  Flüssigkeiten  ist  aus 
Fig.  (3)  ersichtlich.  Auch  die  genaue 
Formel  (74),  wo  c2  nicht  gegen  b ver- 
nachlässigt ist,  liefert  für  Rq  = 10  keine 
merklich  verschiedenen  Werte;  z.  ß.  wird 
nach  (74)  0 = — 0,007  und  nach  (74a) 
0 = — 0,005. 

Da  also  die  Vernachlässigung  von  c2 
gegen  b statthaft  ist,  zeigen  sich  0 und 


Fig.  8. 


|| 

0 

W 

0 

0 

0 

0,5 

— 0,99 

+ 0,05 

1,0 

— 0,94 

+ 0,28 

1,5. 

— 0,74 

+ 0,53 

2,0 

— 0,39 

+ 0,64 

2,5 

— 0,09 

+ 0,58 

3,0 

+ 0,09 

+ 0,43 

3,5 

4-0,18 

+ 0,27 

4,0 

+ 0,19 

+ 0,15 

4,5 

+ 0,16 

+ 0,05 

5,0 

+ 0,13 

— 0,006 

5,5 

+ 0,09 

— 0,04 

6,0 

+ 0,05 

— 0,04 

8,0 

— 0,01 

— 0,02 

10,0 

- 0,005 

+ 0,002 

und  damit  alle  unsere  weitern  Ausdrücke  im  wesent- 
lichen von  der  Elastizität  der  Röhrensubstanz  unab- 
hängig; denn  sie  enthalten  nur  den  Poisson’schen  Coeffizienten 
<9,  nicht  b. 

§ 4.  Bewegung  der  Flüssigkeitsteilchen. 

Ueber  die  Geschwindigkeit  u in  Richtung  der  Röhrenachse 
und  w senkrecht  dazu  in  radialer  Richtung,  geben  die  Gl.  79 
und  83  für  allgemeine  Wellen  und  die  Gl.  96  und  99  für  rein 
radiale  Schwingung  der  Röhrenwand  Aufschluss.  Da  w im 
Verhältnis  m : 1 kleiner  ist  als  u,  so  sind  die  radialen  Geschwin- 
digkeitskomponenten der  Flüssigkeitsteilchen  stets  sehr  klein  im 
Vergleich  zu  den  Longitudinalen.  Die  Werte  der  Amplituden 

von  u und  w,  also  die  Maximalwerte 
dieser  Geschwindigkeiten,  für  allgemeine 
Wellen  und  für  rein  radiale  Bewegung 
der  Röhrenwand,  finden  sich  in  den  nach- 
stehenden Tabellen.  In  der  nebenstehen- 
den Zusammenstellung  sind  die  Ampli- 
tuden von  u für  allgemeine  Wellenbe- 
wegung, für  2 Werte  von  Rq  angegeben 
(je  mit  103  multipliziert),  Fig.  4 zeigt 
die  graphische  Darstellung  derselben. 

Geschwindigkeit  u für  allgemeine  Wellen.  Fig.  4. 


r q 

Geschwindigk.  u 

►o 

II  • 

R q = 10 

0 

5,98 

4,97 

2,0 

5,21 

4,97 

4,0 

0,94 

5,05 

6,0 

— 

5,29 

8,0 

— 

4,96 

10,0 

— 

1,06 

Die  unerwartete  Form  dieser  Kurven,  wie  auch  der  spätem 
Fig.  5 u.  6,  lässt  sich  durch  folgende  Ueberlegung  verstehen. 
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Durch  die  mit  konstanter  Amplitude  erzeugte  Schwingung  der 
elastischen  Umhüllung  werden  die  anliegenden  Flüssigkeitsteilchen 
ebenfalls  in  gleichen  Schwingungszustand  gebracht  und  die  kine- 
tische Energie  wird  nun  auf  die  benachbarten  Flüssigkeitsschichten 
übertragen.  Infolge  der  freien  Beweglichkeit  der  Flüssigkeit, 
können  dieselben  mit  zunehmender  Geschwindigkeit  hin  und  her 
pendeln,  so  dass  bei  dünnen  Röhren,  resp.  kleiner  Frequenz  oder 
grosser  Viskosität  (Rq  = 4)  (gestrichelte  Kurve)  die  Amplitude 
der  Geschwindigkeitskomponente  u bis  zur  Röhrenachse  stetig 
zunimmt.  Es  nimmt  aber  mit  wachsender  Geschwindigkeit  und 
besonders  mit  wachsendem  Geschwindigkeitsgefälle  in  radialer 
Richtung  auch  die  Reibung  zu,  so  dass  diese  Amplituden  bei  ge- 
nügend grossen  Röhren  und  bei  genügend  kleiner  Viskosität, 
wo  u rasch  wachsen  kann,  also  auch  ein  starkes  radiales  Gefälle 
besitzt,  nicht  dauernd  zunehmen,  sondern  ein  Maximum  erreichen 
und  dann  wieder  abnehmen.  Für  Rq  = 10  (ganz  ausgezogene 
Kurve  Fig.  4),  für  eine  Röhre  mit  relativ  grossem  Durchmesser 
resp.  für  grosse  Frequenz  oder  relativ  kleine  Viskosität,  tritt 
diese  Erscheinung  deutlich  zu  Tage.  Im  Abstand  von  ca.  3/5  R 
ist  dieses  Maximum  erreicht,  nachher  nimmt  u bei  Annäherung 
an  die  Röhrenachse  wieder  ab.  Bei  zäherer  Flüssigkeit  dagegen 

zeigt  die  nebenstehende  Tabelle  und  die 
schwach  punktierte  Kurve  der  Fig.  4 für 
Rq  = 4 (welche  für  dieselbe  Röhre,  mit 
dem  gleichen  Radius  R wie  diejenige  für 
Rq  = 10,  nur  mit  veränderlichem  q, 
speziell  für  grössere  Viskosität  gilt),  dass 
in  der  Nähe  der  Röhren  wand  die  Ge- 
schwindigkeitszunahme viel  langsamer  er- 
folgt, also  auch  die  Reibung  viel  geringer 
wird  als  für  Rq  = 10.  Infolgedessen  kann  die  Geschwindig- 
keit mit  Annäherung  an  die  Röhrenachse  immer  noch  zunehmen, 
während  sie  für  die  weniger  zähe  Flüssigkeit  bereits  ihr  Maxi- 
mum erreicht  hat.  Es  ist  also  möglich,  dass,  wie  Fig.  4 zeigt, 
für  zähe  Flüssigkeiten  die  Kurve  bis  zur  Röhrenachse  stetig 
zunimmt. 

Genau  dieselben  Bemerkungen  gelten  für  rein  radiale 
Schwingungen  der  Röhrenwand. 


r 

Geschwindigk.  u 

Rq  = 4 

Rq  = 10 

0 

5,90 

4,97 

0,25 

5,85 

4,97 

0,50 

5,20 

5,14 

0,60 

4,70 

5,29 

0,80 

3,20 

4,96 

1,00 

0,94 

1,06 

Numerische  Angaben  für  allgemeine  Wellenbewegung. 
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Für  die  Amplituden  der  Geschwindigkeit  w,  senkrecht  zur 
Röhrenachse  gelten  analoge  Betrachtungen,  nur  dass  hier  ihr 
Betrag  auf  der  Röhrenachse  Null  ist  und  an  der  Röhrenwand 
der  Maximalwert  erreicht  wird. 

Um  die  Bewegung  der  einzelnen  Flüssigkeits- 
teilchen genau  zu  verfolgen,  genügt  es  nicht,  die  Amplituden 
von  u und  w zu  kennen,  es  spielen  die  auftretenden  Phasen- 
differenzen d u.  d'  gegenüber  der  Schwingung  des  Druckes  eine 
beträchtliche  Rolle.  Noch  deutlicher  ergibt  sich  diese  Bewegung 
durch  Berechnung  der  Elongationen  x und  r der  einzelnen  Teilchen 
für  verschiedene  Zeiten.  In  den  vorstehenden  Tabellen  sind 
für  einen  bestimmten  Fall,  Rq  = 10,  diese  Werte  berechnet, 
wobei  die  Zahlen  mit  103  resp.  104  erweitert  sind.  Dabei  wurde 

m = — 0,005  gesetzt,  was  bei  einer  Frequenz  o — 10,  einer 

Fortpflanzungsgeschwindigkeit  c — 20  m/sek,  entspricht.  Fig.  5 
zeigt  die  betreffenden  Verhältnisse  für  allgemeine  Wellen  für 
die  Mittellage  x = o in  graphischer  Darstellung  für  die  Werte 


2 T 

8 T 

/7T\ 

4 T 

5 T 

6 T 

7 T 

f15T) 

8 T 

8 ’ 

8 ’ 

\ 16/ 

’ 8 ’ 

8 ’ 

8 ’ 

8 ’ 

V 16  J 

’ 8 

(wo  T = Schwingungsdauer)  R — 1,  o — 10  und  v — 0,1. 

In  Fig.  6 sind  die  nämlichen  Verhältnisse  für  rein  radiale 
Schwingung  der  Röhren  wand  aufgetragen. 

§ 5.  Bemerkungen  über  die  Anwendung  der  Formeln 
auf  Bestimmung  derFortpflanzungsgeschwindigkeit 
der  Pulswellen  in  den  Arterien. 

Ueber  die  Fortpflanzungsgeschwindigkeit  der  Puls  wellen 
im  menschlichen  Körper  wurden  seit  Weber  von  Prof.  Münzer19) 
u.  a.  zahlreiche  Versuche  gemacht;  die  Werte  von  Münzer  vari- 
ieren zwischen  5 — 26  m/sek.  je  nach  dem  pathologischen  Zu- 
stand des  Patienten  und  der  Anzahl  Pulsbewegungen.  Um 
unsere  Formeln  auf  die  experimentellen  Resultate  an  wenden  zu 
können,  sollte  in  erster  Linie  der  Elastizitätsmodul  E der 
Arterien  bekannt  sein,  worüber  aber  keine  bestimmten  Angaben 
vorliegen  und  überhaupt  schwer  zu  erhalten  sein  werden,  da  die 
Arterien  nicht  das  einfache  Verhalten  eines  isotropen,  elastischen 
Körpers  zeigen  und  zudem  die  sie  umgebenden  Fettmassen,  Bind- 
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gewebe  etc.  den  Elastizitätsmodul  beeinflussen.  Als  Röhren- 
durchmesser R käme  der  mittlere  Arteriendurchmesser  zwischen 
den  beiden  Messpunkten,  zwischen  denen  die  Fortpflanzungs- 
geschwindigkeit c ermittelt  wird,  in  Betracht,  und  in  entsprechen- 
der Weise  müsste  ein  Mittelwert  für  die  Dicke  h der  Arterien 

bestimmt  werden.  Da  aber  bei  den  Arterien  das  Verhältnis 


nicht  klein  ist,  so  fallen  dieselben  wieder  etwas  ausserhalb  des 
Gültigkeitsgebietes  unserer  Formeln.  Immerhin  ist  es  möglich 

h • E 

in  angenäherter  Weise,  aus  Gl.  (67)  c2  = ^ ^ K,  bei  bekannten 

Dimensionen  der  Arterien  und  bei  bekannter  Viskosität  v einen 
Elastizitätsmodul  E der  Arterien  zu  berechnen.*)  Aber  diese 
Konstante  E hat  natürlich  aus  den  angeführten  Gründen  keine 
einfache  physikalische  Bedeutung  mehr.  Dagegen  zeigen  unsere 
Ausdrücke,  wie  die  Wellengeschwindigkeit  c (wie  auch 
die  Geschwindigkeitskomponenten  u und  w der  einzelnen  Blut- 
teilchen) von  dem  Arteriendurchmesser,  der  Pulsfrequenz 
und  der  Viskosität  des  Blutes  abhängen,  woraus  vielleicht 
Konsequenzen  für  die  Physiologie  gezogen  werden  können.  Ins- 
besondere liegt  auch  die  Möglichkeit  vor,  die  von  Prof.  0.  Frank20) 
ausgeführten  Messungen  an  Manometerschläuchen  auf  ihre  Ab- 
hängigkeit von  der  Zähigkeit  der  Flüssigkeit  hin  zu  prüfen. 

Berechnungen  über  die  Dissipation  von  Energie  infolge  der 
Reibung  wurden  in  der  vorliegenden  Arbeit  nicht  durchgeführt. 
Ueber  direkte  Energie-Messungen  der  Pulswelle  sei  auf  die 
Arbeiten  von  Dr.  Th.  Christen,21)  Prof.  Dr.  Sahli22)  u.  a.  verwiesen. 


*)  Umgekehrt  kann  bei  irgend  einer  elastischen  Röhre,  deren  Dimen- 
sionen und  Elastizitätsmodul  bekannt  sind,  nach  experimenteller  Bestim- 
mung von  c die  Konstante  K berechnet  werden,  aus  welcher  sich  nach 


gegebener  Frequenz  a der  kinematische  Reibungskoeffizient  r 
berechnet  werden. 
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besonders  den  Herren  Prof.  Dr.  J.  H.  Graf,  Prof.  Dr.  G.  Huber 
und  Prof.  Dr.  Förster  für  die  klaren  und  anregenden  Vor- 
lesungen. 
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